Semiflexible Polymer Confined in Closed Spaces 
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We analyze static properties of a strongly confined semiflexible polymer, i.e. either trapped in a 
closed space or compressed by external forces, in an athermal solvent. Like a flexible polymer case, 
we can resort to an analogy with the semidilute solution, but a complication due to the additional 
length scale arising from the chain rigidity results in more diverse behaviours depending on system 
parameters. For each regime, scaling forms of the excess free energy of the confinement are derived. 
Effects of the confinement geometry and the solvent quality are briefly discussed. 
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I. INTRODUCTION 

Macromolecules are very soft objects. As a conse- 
quence, they can be squeezed into spaces with spatial 
size substantially smaller than their preferable size in 
bulk P, H, H, 0]. Understanding behaviours of such 
confined polymers serves as the basis for range of phe- 
nomena encountered in polymer industry, biotechnology 
and various molecular processes in living cells. Thanks 
to the progress in nanofabrication techniques and single 
molecule observations, experiments in this filed are now 
rapidly developing @, @, 0, g, 1, [T(| . 

We have recently pointed out that polymer chains 
can be intensively compressed in certain situations so 
that the excess free energy of the confinement F ex- 
hibits a nonlinear dependence on the chain length N. 
This led us to propose two distinct confinement regimes, 
i.e., weak confinement (WCR) and strong confinement 
(SCR) 11]. The qualitative difference between these two 
are easily understood using the example of a flexible lin- 
ear chain. Suppose that a long flexible chain with the 
bulk size R ~ bN 3 ' 5 (b is the monomer size) is con- 
fined between the parallel plates with their separation 
D < R. In a length scale smaller than D, the be- 
haviours of the chain is not perturbed by the wall, and 
this introduces the blob size D, inside which there are 
g = (D/b) 3 / 5 monomers. The local monomer volume 
fraction inside blob <f> ~ b 3 g/D 3 ~ (b/D) 3 / 4 is indepen- 
dent of the chain length. This corresponds to what we call 
WCR. The extensive nature of the free energy, therefore, 
results in the conventional scaling form of the confine- 
ment free energy F ~ Jb/D) 3 ^ 5 N , that is proportional 
to the chain length 0, Q . The same argument is appli- 
cable for the chain confined in the capillary, too. How- 
ever, it is no longer hold for the chain trapped in the 
closed cavity, where the volume fraction increases with 
the chain length [12]. This provides the simplest exam- 
ple of the SCR, in which several distinct properties are 
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expected [111 ]. 

In ref . [ll[ , we generalized this concept to the branched 
polymers, where the constraint imposed by the connec- 
tivity results in nontrivial confinement behaviours. An- 
other requisite extension is examining the effect of the 
stiffness. The reason is immediately recognized by the 
fact that many biopolymers such as DNA and various 
proteins have a quite large stiffness aspect ratio p = l/b, 
where I is the Kuhn segment length. Even for less stiff 
chains than DNA (note that any real polymers possessed 
some degree of intrinsic stiffness), such a consideration 
may become apparent in denser system, in particular in 
the confined system, where the additional length scale of 
the confinement D would produce a variety of the be- 
haviours. 

It would be convenient to introduce the distinction be- 
tween WCR and SCR for semiflexible chain, too. Scaling 
laws describing properties of a semiflexible linear chain 
trapped in the capillary or the slit (WCR) are the same 
as those of a flexible one if the size of the confinement is 
much larger than the persistence length l p (^ I); the ef- 
fect of the stiffness shows up just through the definition 
of link size. At D < I, however, the chain is obliged to be 
bent sharply. This leads to the introduction of the Odijk 
deflection length lo ~ D 2 / 3 l]/ 3 , which provides / depen- 
dences of various quantities [13]. Recent observation of 
a single DNA chain in a narrow channel have shown the 
crossover from a flexible to a semiflexible chain regime 
around the channel size D ~ I [l(| ■ 

A congeneric level of argument seems to be lacking for 
the semiflexible linear chain in SCR, e.g. in a closed cav- 
ity, and it is the purpose of the present paper to provide 
a qualitative guide for such a problem. In particular, 
we discuss scaling forms of the confinement free energy, 
which is one of the most important quantities in the prob- 
lem, with their physical origin in various regimes. To do 
so, we first discuss the static property of a flexible linear 
polymer in SCR in some details in Sec. [TTJ Based on 
the semidilute solution analogy, we show that the con- 
finement free energy scales F ~ TV 9 / 4 under good solvent 
condition [28| , while F ~ N 3 at "lower" temperature and 
high concentration . Then, in Sec. IIII1 we proceed to the 
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more general case of a semiflexible linear polymer. As we 
shall show, the presence of multiple length scales, i.e., the 
cavity size, the Kuhn length and the correlation length of 
the concentration fluctuation, leads to rich confinement 
regimes even in the athermal solvent. Some of these can 
be viewed as a semidilute solution as in the flexible chain 
case, and in such cases, we can directly apply the result 
known in the theory of semidilute solution of semiflexible 
chains [H, However, there are other regimes, which 
manifest intrinsic effects of the confinement originated 
from either a reduction in the conformational entropy 
or a sharp bending. Throughout the article, we assume 
that the segment concentration of the confined chain is 
not very high, in which the model independent univer- 
sal behaviours can be discussed. To make the discussion 
simple, we also assume that the confinement geometry is 
a sphere (the effect of the shape will be addressed in Sec. 

El). 



II. A FLEXIBLE CHAIN IN SCR 

Let us consider a flexible linear chain (contour length 
L) dissolved in a solvent. Large-scale properties of the 
chain is conveniently modelled as a sequence of N links 
of size b (volume v ~ b 3 ), i.e., the standard bead-on- 
filament model, with the distance between neighboring 
beads along the chain a [l2| . For a flexible chain, the 
natural choice is N = L/b and a = b. The excluded 
volume effect becomes important when the parameter 
z(N) = BN 1 / 2 /a 3 is larger than unity, where B = vt 
is the second virial coefficient. The solvent quality is 
designated by a reduced temperature r = (T — 0)/9 with 
T and 9 being the absolute temperature and a 9 temper- 
ature, respectively. From the condition z(gth) — 1, the 
number of monomers gth = t constituting the ther- 
mal blob is obtained. Inside the thermal blob, the chain 
conformation is Gaussian, leading to the thermal blob 
size £th — br^ 1 . The chain conformation in a bulk solu- 
tion can be envisioned as a self-avoiding walk of thermal 
blobs of size £th — b/r, which leads to the overall chain 
size R ~ WV 3 / 5 !- 1 / 5 . 

Then, consider that the chain is brought into a spher- 
ical cavity of size D < R. This corresponds to the sim- 
plest example of the SCR. A key observation is that the 
strongly confined chain can be viewed as a semidilute so- 
lution with the volume fraction <p ~ b 3 N/D 3 [ll[, i.e., the 
correlation of concentration fluctuations is suppressed at 
the length scale 
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The above expression for the correlation length is ob- 
tained from a scaling argument @ (by imposing £ must 
be (i) independent of N and (ii) equal to R at the overlap 
concentration) and valid as long as 4> < T £ > £th- 

This corresponds to the fluctuating semidilute regime. In- 
side the blob, the correlation is evident, thus, the number 



of monomers in each blob is obtained from the relation 
£ ~ &5 3 / 5 t 1 / 5 . At larger length scale, the application of 
the mean-field argument is valid with the blobs of size £ 
being the renormalized monomers. A free energy arising 
from volume interactions is , thus, evaluated as 



F, 



vol 



AM (2) 



where B* ~ £ 3 is the renormalized second virial coeffi- 
cient. This is equivalent to assigning the energy on the 
order of fc^T to each blob. 

In the opposite case, 4> > t ^=^> £ < £ t h, the 
chain conformation is Gaussian in all the length scale, 
which indicates the applicability of the mean field the- 
ory. This corresponds to the semidilute in 9-solvent 
regime. The volume interaction energy per unit vol- 
ume can be expressed in a virial expansion f vo i/kBT ~ 
l/2Bc 2 + l/3Cc 2 +, where B = vt, C = v 2 are the second 
and third virial coefficients and c = <f>/v is the monomer 
concentration. The monomer pair-correlation function is 
calculated by a random phase approximation [H, 0, [H[ 
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where, the correlation length in this regime is, instead of 
eq. ((J) , given by 
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^ MF ~ 12 \k B T dc 
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where the osmotic pressure (due to volume interactions) 
can be expanded in a series of power c as H/(ksT) = 



Bc 2 + 2Cc 



Under the present condition <fi > r, the 



interaction is dominated by the three body contribution, 
therefore, the total volume interaction contribution to 
the confinement free energy and the correlation length 
are evaluated as 



k^f 
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(5) 
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Note that short chains (N < g t h) are in the 9 solvent 
regime already in the bulk isolated state, to which eq. 
(J5]) and $Q always hold upon the confinement. The con- 
finement regimes of a flexible chain in a closed cavity is 
presented in diagrammatic form in Fig. [TJ 

The scaling form of eq. ([2]) in the athermal solvent 
has been confirmed by a recent computer simulation [1 61 ] , 
while the conjecture has been proposed [ll| that the form 
of eq. ([5]) would be relevant to resolve controversial exper- 
imental results on the partitioning coefficient of a flexible 
polymer into a protein pore H [13, [H, EH ■ 

In both regimes, aside from the volume interaction 
free energy (bulk term), there is a contribution from the 
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FIG. 1: A diagram of confinement regimes for a flexible chain 
in a closed sphere (logarithmic scale) . Region (I) and (II) cor- 
respond to fluctuating semidilute and semidilute in 8-solvent 
regime, respectively. Shaded areas are irrelevant, i.e., upper 
left is a bulk (D > R) and lower right is not accessible (0 > 1). 
Border lines between regimes are (a) D/b ~ r 1/5 N 3/B , (b) 
D/b ~ (A7r) 1/3 , (c) D/6 ~ A^ 1/3 , (d) D/b ~ Af 1/2 . 



nonuniform link concentration AF (surface term) [Tl| 
These are calculated in Appendix [5] and given by 



AF 
~kBT 



— T*N* 



in the fluctuating semidilute regime and 



AF 



4>D 



N 



(8) 



in the semidilute in ^-solvent regime. 

For a flexible chain, the dominant term is always given 
by the volume interaction, i.e., F ~ F vo i. However, as 
we shall see below, there is a regime, in which AF plays 
a central role for a semiflexible chain. 



III. A SEMIFLEXIBLE CHAIN IN SCR 

Our result is summarized in the diagram of Fig. [2] To 
discuss this diagram, we first assume (i) D > I and (ii) 
L > pD (here p is the stiffness ratio mentioned in Sec. 
|T|. As we shall see later, these conditions are requisite 
for a confined chain to be viewed as an analogue of bulk 
semidilute solution. 

Let us begin with the evaluation of the volume inter- 
action parameter z = BN 1 / 2 /a 3 . There is a freedom 
on how to break the chain into links, and we adopt the 
simplest, in which the chain is represented as N = L/b 
links with their volume v ~ b 3 (6 is the chain thick- 
ness). Then, for the model to be consistent, the dis- 
tance between links should be set as a = (bl) 1 / 2 [r3 |. 
It should be emphasized that results obtained are inde- 
pendent of such a modelling (one can also set N = L/l, 



FIG. 2: A diagram of confinement regimes for a semiflexible 
chain in a closed sphere under athrmal condition (logarithmic 
scale). Each labelled region designates (I) fluctuating semidi- 
lute, (II) mean-field semidilute, (III) liquid crystalline, (IV) 
ideal chain and (V) bending regime, respectively (see the main 
text for details). Shaded areas are irrelevant, i.e., upper left 
is a bulk (D > ii) and lower right is not accessible (<j> > 1). 
Border lines between regimes are (a) D/b ~ p 1/5 N' i/5 , (b) 



D/b ~ pN 1/3 , (c) D/b 
D/b ~ (piV) 1/2 , (f) D/b z 



- (pN) 1/3 , (d) D/b ~ Ar 1/3 , (e) 
N/p, (g) D/b ~ p, (h) D/b ~ N. 



v ~ bl 2 and a = I). For simplicity, we shall focus on 
the purely repulsive case, i.e., athermal solution: r = 1 
(thus, B = v). From z = 1, we find the length scale 
£th — bp 2 , below which the excluded volume effect is 
insignificant (the corresponding number of monomers is 

1/2 

ag^ ~ £,th)'- for semiflexible chains, £th ^ b even in the 
athermal limit. As a consequence, the chain with N < p 3 
behaves as a Gaussian chain, i.e., its size in bulk is given 
by R ~ {IL) 1 / 2 . On the other hand, the excluded volume 
interaction should be taken into account for the longer 
chain, i.e. R ~ & h (N / 'g th ) 3 ^ ~ {bl)^ 5 L 3 / 5 . 

Along with the same argument as in the flexible chain, 
one can distinguish (I) fluctuating semidilute regime and 
(II) mean-field semidilute regime depending on the rel- 
ative ratio of t; t h and the correlation length of density 
fluctuations. 

(I) Fluctuating semidilute regime: pN 1 ^ 3 < D/b < 

(p)l/8jV8/S 

According to the analogy with the semidilute solution, 
a confined chain is viewed as a dense packing of blobs of 
size £ = £,th(g/ gth) 3 ^ 5 (g is the number of monomers in 
each blob). This indicates the condition g/£ 3 = N/D 3 , 
from which one finds 



p * N 4 



(9) 



The dominant term in the confinement free energy comes 
from the volume interaction 
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(II) Mean-field semidilute regime: (pN) 1 ^ 3 < D/b < 
minlpiV 1 / 3 ,^ 1 ^} 

A correlation length £ becomes shorter with the in- 
crease in L (or decrease in D) and becomes comparable 
to £ th at D ~ L/V 1 /^ ~ p -3)_ xhis indicates a 

crossover to the mean-field semidilute regime. One can 
also recognize the meaning of this condition through the 
Ginzburg criterion: 



d 3 rd 3 r'(Sc{r)Sc{r')) < c 2 



^rj_my < 

\k B Tdc) ~ 

^P~ 3 < 



cvp 2 



(11) 

(12) 
(13) 



where bracket denotes the statistical average and ~ 
£mf 1S t ne correlation volume. The correlation length 
is calculated from eq. ((4|) with ignoring the third and 
higher order terms in virial expansion 



MF 



(14) 



(Note that the eq. I|12p also provides the borderline for a 
flexible chain discussed in the previous section.) 

Here again, the confinement free energy is dominated 
by the volume interaction, which is evaluated by the 
mean field theory. In the case under consideration 
((ML) 1 / 3 < D), the binary interaction is the most dom- 
inant than higher order terms in the virial expansion, 
thus, 



F 



kuT knT 
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This scaling form has been proposed on the basis of a 
self-consistent field calculation [20], but here one should 
notice the range for the applicability, i.e. only valid for 
a semiflexible chain in the present regime. 

(Ill) Liquid crystalline regime: p < D/b < (pN) 1 / 3 
If the density is further increased, the volume inter- 
action becomes higher and the correlation length be- 
comes comparable to the Kuhn segment length at D ~ 
(&/L)V3( ^ ~ p -iy T hi s indicates some role of 

model specificity on the flexibility mechanism at higher 
concentrations, analysis of which is beyond the scope of 
the present discussion. Here, we just note what is ex- 
pected from the bulk theory. According to the standard 
theory [2l| , the system responds by breaking the isotropic 
symmetry, i.e., the competition between steric repulsion 
(evaluated by the second virial approximation) and orien- 
tational entropy results in the first-order phase transition 
at the volume fraction <fi ~ p~~ x . The transition proceeds 
with the coexistence of isotropic and nematic phases at 
c\p~ x < <j) < C2p~ 1 (numerical coefficients depend on 
the flexibility mechanism), which suggests an interesting 
possibility of an "intra-chain" separation for a confined 
chain. The free energy in the nematic state is dictated by 



the loss of the orientational entropy, thus, evaluated to 
be proportional to the number of the statistical segment 
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This is just a crude approximation since the degree of the 
ordering increases (even rather slowly) with the concen- 
tration. At much higher concentration, this would not 
be certainly correct due to higher order terms. 

(IV) Ideal chain regime: max{p,p _1 A^} < D/b < 
tpN) 1 ' 2 

We now consider the meaning of the second condition 
p~ l L < D necessary for the bulk semidilute solution 
analogy. To do so, let us compare the correlation length 
of the concentration fluctuation with the cavity size. In 
the fluctuating semidilute regime, the condition D > £ 
(eq. [5]) is equivalent to trivial condition for the confine- 
ment D < R. On the other hand, in the mean-field 
semidilute regime, D > £,mf (eq. 1 14[) leads to nontrivial 
condition D < p~ x L. Note that for chains not too long 
(N < p 3 ), this length p~ x L is still smaller than the nat- 
ural chain size in bulk R = {IL) 1 / 2 . Here the correlation 
length calculated in bulk exceeds the system size, indi- 
cating that the confined chain is no longer analogous to 
the bulk semidilute solution at D > p~ x L. This regime 
can be regarded as a critical point analogue in a finite 
confined system ("critical region") and never expected 
for a flexible chain. Approaching to the line D = p~ l L 
from the mean-field semidilute regime, the bulk region 
decreases and finally the whole system becomes nonuni- 
form. The free energy associated with the nonuniformity 
is no longer regarded as the surface term and the confine- 
ment free energy is evaluated from the integral of eq. (|A4[) 
by replacing £ with D. 

( A\ and A2 are numerical coefficients of order unity) . In 
the present case with D > p L, the first term (aris- 
ing from the reduction in the conformational entropy) 
is dominant, i.e., excluded volume interactions are dis- 
regarded and the confinement free energy is approxi- 
mated by that of an ideal chain trapped in a cavitiy 
F/(k B T)~(b/D) 2 pN. 

(V) Bending regime: D/b < min{p, N} 

The behaviours of the confined chain should be 
strongly dependent on the flexibility mechanism once 
the cavity size becomes smaller than the Kuhn segment 
length. In particular, it is obvious that a freely-jointed 
chain cannot be placed inside such a small space. We 
consider here a worm-like-chain, which is suitable as a 
model of DNA packing inside virus capsid. A worm-likc- 
chain possesses a uniform elastic modulus along the chain 
contour k = ksTl/2. When such a chain is confined in a 
small cavity with size D < I, each chain section is forced 
to retain the curvature on the order of D . In average, 
the bending energy is almost uniformly distributed along 
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the chain. Thus, our estimate for the confinement free 
energy is 



Fbe 



L 



F 



k B T k B T k B T D 2 



pN 



(18) 



The loading force is not length dependent, which is sup- 
ported by the result of computer simulation if not the 
segment concentration is too high (<j> < 0.2) [22|]. Again, 
at high concentrations, segmental interactions become 
important so that the confined chain may exhibit some 
structural order and eq. (|18p is no longer applicable. 
There are several studies for such a situation in connec- 
tion to the problem of DNA packing inside virus cap- 
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IV. DISCUSSION 

In this section, we shall discuss some generalization of 
the obtained results. We have assumed the confinement 
in a spherical container only. This is just for the sim- 
plicity and the generalization to asymmetric geometry 
is straightforward. For instance, let us compare a chain 
confined in a rectangular box of volume V — L\LiLj, 
[L\ < L 2 < L 3 ) and that in a sphere with equal volume 
(radius D = (3V r /47r) 1/ ' 3 ). Since the link concentrations 
are the same, the confinement free energies are also the 
same in the leading order in regime (I) and (II), and the 
slight difference just comes from the surface term AF. 
We expect that this explains the subtle shape effect of the 
capsid found in a recent simulation that a sphere packs 
more quickly and ejects more slowly a flexible chain than 
an ellipsoid [25|. On the crossover from the confined to 
the unconfined situation, there appears a WCR between 
SCR (£ < Li) and bulk (R < Li), where the correlation 
length is set by the smallest size (£ = Li), i.e., a chain 
confined in a slit. The border between regime (II) and 
(IV) is given by LiL^fLx ~ L/p, and in regime (IV) the 
entropic contribution as a dominant term is set by the 
smallest size, i.e., F/(k B T) ~ (b/L^pN. 

Effect of the solvent quality is also easily included by 
adding the reduced temperature term in the second virial 
coefficient, i.e., B — vt, provided that the corresponding 
scmidilute solution stays in the one-phase region. How- 
ever, for a polymer with large p, the (intra-chain) phase 
separation and/or the transition to highly ordered ne- 
matic states would occur above 9 point as expected from 
the theory of the semiflexible polymer solution [12, Hl| . 



V. SUMMARY 



and the cavity size. Three of them have semidilute so- 
lution analogues, but other two are intrinsic to a con- 
fined chain. We have not aimed to describe high concen- 
tration states (liquid crystalline and bending regimes), 
since properties there are not universal and more elabo- 
rate models are required [12, HU [H[ . For other regimes, 
the obtained results are rather general thanks to the uni- 
versality of non-concentrated polymer solutions. 

Approaching to the flexible chain limit (p — > 1), all the 
regimes merge into a fluctuating semidilute regime, in 
which the confinement free energy is given by F ~ N 9 ^ 4 . 
Flexible chains such as polyethylene glycol belong to this. 
Note that for a confined chain, the presence of a slight 
segmental attraction leads to a semidilute in 0-solvent 
regime (F ~ N 3 ). 

Now we discuss the experimental accessibility of the 
present predictions. As an example of semiflexible chains, 
let us take up DNA molecules. The Kuhn length and the 
width of the double strand are I ~ 100 nm and b ~ 2 
nm, respectively, in usual aqueous conditions, thus, the 
aspect ratio p ~ 50 (note that the electrostatic repulsion 
results in some thickening at low ionic strength). First, 
we see that a fluctuating semidilute regime requires very 
long DNA [L > bp 3 ~ 250^m). Although many nat- 
ural DNA in biological origin is long enough to satisfy 
this, it is difficult to treat such long chains in a reliable 
manner. Therefore, other regimes are more easily acces- 
sible and relevant to most experiments. For moderate 
(but still rather long) chains (5/im < L < 250/im), a 
variety of regimes are expected depending on the cavity 
size. For example, if T4 DNA (L ~ 56/zm) is trapped 
in the cavity with size ca. 1/Ltm, a mean-field semidilute 
regime (F ~ N 2 ) is most probable, while an orientation- 
ally ordered state may be observed in the small cavity of 
a few submicron size. In the case of A DNA (L ~ 16/im), 
an ideal chain regime (F ~ N) is expected as long as 
the cavity size exceeds a few submicron, etc. For shorter 
chains (L < 5/im), we expect an ideal chain regime, or a 
bending regime depending on the relative size ratio be- 
tween I and D. 

Although qualitative (all the numerical coefficients of 
order unity are not determined), we believe that the 
present discussion provides a useful guide for emerging 
techniques operating long polymers, such as DNA, in 
nm ~ am scale spaces, where molecules may be highly 
compressed geometrically 0] or by external field [1]. It 
is also relevant to phase behaviours of composite soft 
matter systems such as a mucroemulsion + polymer sys- 
tem (26j . In certain situations (e.g., a polymer com- 
pressed by electric filed against a wall), a nonuniform 
link distribution should be taken into account [2, [27| . 



We have discussed static properties of a long linear 
polymer confined in the closed cavity. Guided by an anal- 
ogy between a long confined chain in SCR and semidilute 
solutions, we have found five distinct regimes for a semi- 
flexible chain (p ^> 1) depending on the contour length 
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APPENDIX A: CALCULATION OF AF 



for the depletion layer near the wall. The essentially same 
form as eq. IA3I is given in ref. [|[ for the monomer con- 
centration profile of semidilute solution near the repulsive 
wall. The surface free energy is evaluated as 



Here we review briefly what is needed for our discus- 
sion. Consider a long semiflexible chain confined in a 
spherical cavity (radius D). The solvent is assumed to 
be athermal and the cavity wall is repulsive (no adsorp- 
tion). The local link concentration is almost uniform 
(c(r) = c ~ N/D 3 ) within the cavity, but drops the near 
the wall toward c(r = D) = 0. If the local fluctuation is 
negligible, the concentration profile can be obtained by 
self-consistent field calculation Q. The self-consistent 
equation is 
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Uscf{t) 
k B T 



ip(r) = eijj{r) 



(Al) 



where ip(r) is the grand state eigen function of the prop- 
agator (normalized as c(r) = |-0(r)| 2 ) and the potential 
is 



v\ip(r)\ 



(A2) 



By imposing the boundary condition ip(D/2) = and 
ijj{r) = c 1 / 2 far away from the wall, one finds the grand 
state eigen value e = vc (consistent with the free energy 
eq. (US])) and the profile 



c(r) = c tanh 2 



D 



(r < D) 



(A3) 



where £ = b{p/'&(j)) 1 t 2 , which is, in essence, the correla- 
tion length (( ~ £mf from eq. Q3])), signifies the range 



AF 

k^f 



D 

dS I dr 

D-C 



(> 



(W>M) 2 + ^#MI 4 K A4 ) 



(A5) 



bp 2 D< 
T 3 



The same result is obtained from a slightly different way. 
The surface free energy is equal to the work required to 
create a depletion layer against the osmotic pressure. In 
a mean-field approximation, the osmotic pressure is given 
by 



n = fc s i% 

b 6 



(A6) 



Therefore, AF ~ j dS (II x £mf), which is shown to be 
equivalent to eq. (|A5|) . 

The latter approach can be conveniently applied to 
other situations such as a 9 solvent and a fluctuation 
dominant regime. In the semidilute in 0-solvent regime, 
the osmotic pressure is approximated as II ~ kBT(f> 3 /b s 
(dominated by third virial term). The surface free en- 
ergy for a flexible chain (p = 1) is, thus, AF / {ksT) ~ 
(D(j)/b) 2 ~ (D/^mf) 2 (eq. ©). In the fluctuating 
semidilute regime, II ~ ^T/^ 3 , where one should notice 
the correlation length is obtained from the scaling argu- 
ment (eq. O,©). This leads to AF/(k B T) ~ (D/£,) 2 = 

^3/2 (j3r) l/2 ( £ /5) 2 
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